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ABSTRACT 

Elliptical galaxies are modelled as Sersic luminosity distributions with density profiles 
(DPs) for the total mass adopted from the DPs of haloes within dissipationless ACDM N- 
body simulations. Ellipticals turn out to be inconsistent with cuspy low-concentration NFW 
models representing the total mass, nor are they consistent with a steeper —1.5 inner slope, 
nor with the shallower models proposed by Navarro et al. (2004), nor with NFW models 10 
times more concentrated than predicted, as deduced from several X-ray observations: the mass 
models, extrapolated inwards, lead to local mass-to-light ratios that are smaller than the stellar 
value inside an effective radius (R e ), and to central aperture velocity dispersions that are much 
smaller than observed. This conclusion remains true as long as there is no sharp steepening 
(slope < —2) of the dark matter (DM) DPs just inside 0.01 virial radii. 

The too low total mass and velocity dispersion produced within R e by an NFW-like total 
mass profile suggests that the stellar component should dominate the DM one out to at least 
R e . It should then be difficult to kinematically constrain the inner slope of the dark matter 
DP of ellipticals. The high concentration parameters deduced from X-ray observations appear 
to be a consequence of fitting an NFW model to the total mass DP made up of a stellar 
component that dominates inside and a DM component that dominates outwards. 

An appendix gives the virial mass dependence of the concentration parameter, central 
density, and total mass of the Navarro et al. model. In a 2nd appendix are given single integral 
expressions for the velocity dispersions averaged along the line-of-sight, in circular apertures 
and in thin slits, for general luminosity density and mass distributions, with isotropic orbits. 
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1 INTRODUCTION 

There has been much recent progress on the photometric character- 
ization of elliptical galaxies. Whereas a variety of models for the 
optical surface brightness profiles of ellipticals h ave been used in 
the past, such as the Hubble-Reynolds law <Revnoldsll9H . the an- 
alytical King 1 1962) or modified Hubble law, the projection of the 
isothermal spheres truncated in phase space jKindll966l) and the 
R 1 / 4 law Ide Vaucouleursl [r948). none provides an adequate fit to 
the surface photometry of the large majority of elliptical galaxies. 
However, there has been a recent consensus on the a pplicability 
to vir t ually all elliptical galaxies jC aon 1 Ca paccioli. & D'Onofrid 
Il993t lBertin.Ciotti.&DelPrincipell2002h of the generalization 
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(hereafter, Sersic law) of the _R 1//4 law proposed bv lSersid fl968). 
which can be written as 



I{R) = Io exp 



(R_\ 1/r 
\asJ 



(1) 



where / is the surface brightness, as the Sersic scale parame- 
ter, and m the Sersic shape parameter, with m = 4 recovering 
the law. Moreover, strong correlations have been reported 

between the shape parameter m and either luminosity or effec- 
tive (half projected l ight) radius R e jCaonetaljfl993l see also 
IPrugniel & Simierl 19971 and references therein). 

On the other hand, there is still much uncertainty on the 
importance of dark matter in elliptical galaxies, especially in 
their outer regions. Even in the spherical case considered in this 
paper, kinematical modelling of ellipticals usually cannot dis- 
entangle the degeneracy between the uncertainty of their grav- 
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itational potential and that of their internal kinematics (the 
anisotropy of their velocity ellipsoid), unless veloc ity profiles or 
at least 4th order ve locity moments are con sidered (Merritt 1987; 

3ll992t Ivan der Marel & Frarudll993l lGerhardll993t 



Rix&WhitL . _ 

Lokas & Mamon 2003; Katgert et al. 2004) 



Analyses of diffuse X-ray emission in elliptical galaxies have 
the advantage that the equation of hydrostatic equilibrium has no 
anisotropy term within it, so in spherical symmetry, one easily de- 
rives the total mass distribution through (e.g. Fabricant et al. 1980) 



M(r) = 



kT(r) 



' dlnre dlnT 
- d In r d In r 



(2) 



Gfirrip 

where k is the Boltzmann constant, while T, n and pm p are respec- 
tively the temperature, electron density and mean particle mass of 
the plasma. However, it is crucial to measure T(r) and its gradi- 
ent, and unfortunately, even with the two new generation X-ray 
telescopes XMM-Newton and Chandra, it is difficult to achieve 
such measurements beyon d some fract i on of 1/2 the v irial radius 
for galaxy clusters I Arnaud et al. 2002; Pratt & Arnaud 2002J) but 
much less for elliptical galaxies. Moreover, the X-ray emission 
from elliptical galaxies is the combination of two components: dif- 
fuse hot gas swimming in the gravitational potential as well as di- 
rect emission from indi vidual stars, and it can be h ighly difficult to 
disentangle the two (see Brown & Bregman 2001). 

On the theoretical front, cosmological simulations of large 
chunks of the Universe dominated by cold dark matter (CDM) 
have recently reached enough mass and spatial resolution that 
there appears to be a convergence on the structure and in- 
ternal kinematics of the bound structures, usually referred to 
as haloes, in the simulations. In particular, the density pro- 
files appear to converge to one with an outer slope of ~ —3 
and an inner slop e between —1 iNava goJrgik J _&_Whitejll^9^, 
19961) and -3/2 |Fukushige & Makinclll997l; iMoore et alJll999t 
Ghigna et al. 2000), see also Power et al. (2003); Fukushige et al. 
120041) . 

In this paper, we consider the general formula that ljing & Sutd 
1 2000) found to provide a good fit to simulated haloes: 



p(r) 



(-r M-) 

V a h J L \a h ) 



(3) 



where a = 1 (hereafter 'NFW') or 3/2 (hereafter 'JS-1.5'), and 'ft' 
stands for halo. These profiles fit well the density profiles of cos- 
mological simulations out to the virial radius r v , wherein the mean 
density is A « 200 times the critical density of the Universe. 1 
The ratio of virial radius to scale radius is called the concentration 
parameter: 



a h 



(4) 



Very recently, a number of studi es have pro posed better ana- 
lytic fits to the density l Navarro et al. 2004 ; Diemand et al. 2004 ) or 
circular velocity iSto ehr et alj2002tlStoehj2 005) profiles of simu- 
lated haloes. In particular, the formula of lNavarro et alj 



din p 
d In r 



r-2 



i/n 



(5) 



: 0.3, Qx =0-7, the approx- 



1 For the standard ACDM parameters, f!„ 

imate formulae of Kitayama &^utd Il99d) and|Br^aii^^jorrnan| Jl998|) 
yield A ~ 100 (see also lEke et al]ll99fl and lLokas & Hoffman! l200ll) 
where the exact value is 101.9. However, many cosmologists prefer to work 
with the value of 200, which is close to the value of 178 originally derived 
for the Einstein de Sitter Universe (f! m = 1, £l\ = 0). 



where p ~ 6 and r_2 is the radius of slope —2, is attractive be- 
cause it converges to a finite central density at very small scales and 
steepens progressively at outer radii to produce a finite mass. More- 
over. lNavarro et al J ob tained their fits t o the logarithmic slope of the 
density profile, while loiemand et"al] fit to the density profile and 
Stoehr and collaborators to the circular velocity profile. Since the 
density and circular velocity respectively involve single and double 
integrals of the logarithmic slope profile, the latter two can mask 
subtle variations only pic ked up in the l ogarithmic slope. For these 
reasons, we shall add the Navarro et al. model in the present work. 

In a previous paper l Lokas & Mamon 2001, hereafter paper 
I), we showed that projected NFW density profiles could, in princi- 
ple, be fit to Sersic profiles with m ~ 3. However, there are three 
reasons to disregard this match as a proof that mass follows light 
throughout elliptical galaxies: 

(i) As shown in Fig. 12 of paper I, projected NFW density pro- 
files fit the Sersic profiles only within a narrow range of Sersic 
shape parameters 2.7 ^ m ^ 4.0, given reasonable NFW con- 
centration parameters with (5 ^ c ^ 22), whereas ellipticals are fit 
by Sersic profiles in the much wider ra nge; from m = 0.5 — 0.6 for 
cluster dwarf ellipticals iCaonet all 199311 Binggeli & Jerie 
iMarauez et alfcOOd) to m = 16 ICaon et all) or 8 iGraharrl 19981) . 
7 ( D'Ono friol200lh or at least 5.6 jMarquez et alj200o¥ 

(ii) The fits produced enormous effective radii, much larger than 
observed. 

(iii) The fits produced m increasing with concentration param- 
eter c (again Fig. 12 of paper I — which plots 1/m versus c). 
But, c is known to decrease w i th mass within the virial radius 
iNavarro. Frenk. & Whitelll997t Ijing & Sutd bOOCi). while m is 
know n to increase with luminosity TcaorTet^l IPrugnTel & Simier] 
1 19971 see Fig.0below). Hence, one arrives at the very unlikely re- 
sult that galaxy luminosity decreases with increasing mass within 
the virial radius. 

The conclu sions of paper I have been confirmed in a recent analy- 
sis of iMerritt et alj 120051) . who have recently shown that the ha- 
los in dissipationless cosmological TV-body simulations are also 
well fit in projection by a Sersic model with m ~ 3. This is 
not surprising since the NFW model provides an adequate fit to 
the density profiles of halos in dissipationless cosmological sim- 
ulations. But, in contrast to the case of the divergent-mass NFW 
models, point |(ii)| should not be relevant for convergent-mass ha- 



los described by the Sersic or lNavarro et alj models. Nevertheless, 
the relevance of points |(T)| and |(iii)l remain relative to the work of 
IMerritt et ail 

Now there are several strong indications that mass does not 
follow light in elliptical galaxies: first, from the kinematics of neu- 
tral gas lBert olaetaiIT9 931. and second from X-rays (eq. El), al- 
beit with the simplifying and probably optimistic assumption of 
isothermali ty. Some studi es point to M/L opt i ncreasing out wards, 
such as l Jones et all 1 19971) for NGC 1399. Also. lBuote et alJ t2002) 
use the variation of the ellipticity of the X-ray isophotes to conclude 
that mass does not follow light in an elongated elliptical galaxy 
(NGC 720). 

Recent ana lyses of X-ray data by ISato et alj j2000), 
Wu & Xudi200bl). andlLlovd-Davies & Ponmad <2002l) (quoted in 



Llwd^avie^t^l]l2o62l) conclude (using eg. 0]) that the total 



mass of ellipticals, groups and clusters are indeed well fit by an 
NFW model, but with typically ten times larger concentration than 
measured in cosmological simulations for the mass range of ellip- 
tical galaxies. Surprisingly, there have been no detailed studies of 
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these high concentration total density profiles arising from the X- 
ray observations. 

Whereas the NFW profile was first established with cos- 
mological simulations including a dissipative gaseous component 
iNavarro et alj Il995l) . most confirmations of the NFW or JS- 
1.5 profiles have come from simulations without gas, which can 
achieve a much greater mass resolution. Hence, one may ask if the 
density profiles of haloes found in cosmological iV-body simula- 
tions of cold dark matter in a flat universe with Q.a — 0.7 (globally 
referred hereafter as ACDM) represent the total density of observed 
cosmic structures or only their dark matter component. Given that 
the average baryon fraction in the Universe is small, one could think 
that when structures such as elliptical galaxies form by collapse, 
the baryons simply follow the dark matter without affecting it, and 
the total mass profile would resemble the low concentration cuspy 
models found in the cosmological simulations. Alternatively, since 
gas dynamics is dissipative, there may be a density threshold be- 
yond which the gas will decouple from the dark matter, and then 
the central regions of ellipticals would be dominated by gas and 
later stars that form from it. The low-concentration cuspy density 
profiles found in the cosmological simulations would then only ap- 
ply to the dark matter. Moreover, if the baryons dominate the dark 
matter at small radii, the dark matter will re-adjust itself within the 
gravitational potential dominated by the baryo ns (in a pr ocess often 
referred to as adiabatic contraction, see, e.g., iGnedin et al]|2004l) . 
so that even the dark matter density profile may differ substan- 
tially from the predictions of the dissipationless ACDM simula- 
tions. Also, if ellipti cals form by major mergers of spirals iToomrd 
Il977t l Mamon 1992), the distribution of baryons (mainly stars) will 
be set by violent relaxation operating during the merger, and if the 
baryonic fraction is low in the outer regions of the two merging spi- 
rals, the merger remnant should also show a lower baryon fraction 
in the outer regions. 

In the present paper, we abandon the hypothesis of paper I 
of constant mass-to-light ratio, and we consider the next simplest 
approach, asking ourselves whether the radial profiles of density 
coming out of dissipationless cosmological TV-body simulations 
are consistent with the observations (surface photometry and spec- 
troscopy) of elliptical galaxies, for the total mass (Sec. 0, and 
whether low or high concentration parameters are required. We will 
not consider here the response of the dark matter component to the 
dissipative baryonic component (e.g., the adiabatic contraction of 
the dark matter). 

We begin, in Sec. [2] with a summary of the luminosity and 
mass models of elliptical gal axies that we a dopt in this paper. In 
a companion paper (M anion & Lokas 2005), we go one step fur- 
ther and confront the general observed trends in elliptical galaxies 
with the predictions from a 4-component model of ellipticals with 
stars, dark matter, hot gas and a central black hole, allowing for 
slight radial velocity anisotropy, as seen in cosmological A-body 
simulations. 



2 BASIC EQUATIONS 

2.1 Distribution of optical light 

The Sersic (eq. Q) optical surface brightness profile that represents 
the projected stellar distribution, c an be deprojected accordin g to 
the approximation first proposed by Prugniel & SimieiJ 1 19971) 

£(r) = hl(r/as) (6) 
J(x) ~ x~ p exp (~x 1/m ) , (7) 



T(2m) 



Jo 

1a s ' 



T[(3-p)m] 
1.0 - 0.6097/m + 0.05463/m 2 



(8) 
(9) 



wher e T(a) is the ga mma function. The last equation is from 
llima Neto et alj \ 19991) who argue that equations J6j and jTJ then 
provide a better deprojection of the Sersic profile (eq. Q), and in 
a wider range of m (good to better than 5% accuracy for 0.55 < 
m < 10 within 0.01 < R/R e < 1000, where over 99.5% of 
the light lies) than a previous approximation of p(m) proposed by 
IPrugniel & Simienl 

The integrated l uminosity c orresponding to equations J6j, 
and {8j is then iLima Neto et alj) 



L 3 (r) = LL 3 (r/a s ) 

7 [(3 - p)m, x l/r - 



L 3 (x) 



r[(3-p)m] 



(10) 

(11) 



where 7(0,3;) is the incomplete gamma function and where the 
total luminosity of the galaxy is 

L = 2%m r(2m) Io a% = AnrnF [(3 — p) m] £1 o| , (12) 

as obtained by lYoung & Currid 1 19941) from the Sersic surface 
brightness profile of equation 0, and whi ch matches exact ly the 
total luminosity obtained by integration of iLima Neto et all 's ap- 
proximate deprojected profile. 

Fig. Q shows how the Sersic parameters are correlated. The 
data are taken from 



(i) iBinggeli & Jerienl ll998l) . who computed Sersic parameters 
from fitting the cumulative luminosity profile for dwarf ellipticals 
in Vir go. 

(ii) iMarauez et alj fcoOO). who did the same for normal ellipti- 
cals in the Coma, Abell 85 and Abell 496 clusters. We converted 
their angular sizes to distances, by simply assigning Hubble dis- 
tances to each cluster, assuming no peculiar motion relative to 
us. We also converted their ^-magnitudes to 73-band luminosi- 
ties, assuming B — V — 0.96 (typical of elliptical galaxies, e.g. 
Fukug ita et all 19951) and (e.g. Colina et al. 1996) 



Mr 



5.45 



(13) 



(iii) iD'Onofriol 1200 ll) . who performed 2D fits of Sersic + ex- 
ponential models for galaxies in the Virgo and Fornax clusters. We 
corrected their luminosities by replacing their unique distance mod- 
ul us by the surface br ightness fluctuation (SBF) distance modulus 
of iTonrv et alj 1200 ll) when available, or otherwise by adding the 
median difference (0.18) between the SBF distance moduli and the 
unique value adopted bv ID'Onofriol for their 23 galaxies in com- 



In principle, the 2D fitting method o flD'Onofrid should be the 
m ost reliable, but for given clusters, the means trends in the data 
of Ma rauez et alj show much less scatter for individual clusters. In- 
terestingly, the mean trends for Coma and Abell 496 agree well 
with the mean trend for the ellipticals and lenticulars o fD' OnofrioL 
while the galaxies in the cluster Abell 85 (open triangles in Fig.0 
are systematically offset from the galaxies of the other two clus- 
ters, in such a manner that its distance, as simply derived from its 
redshift, appears 80% too large (as is clear from the correlations 
of m vs. L, and consistent with the correlations in the other plots 
of Fig. 0, but it is highly improbable that the distance to Abell 
85 is overestimated b y such a large factor . Also, the dwarf ellipti- 
cals analysed bv lBinggeli & Jerienl ^98) appear to have a shape- 
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10 9 10 10 10 11 10 9 10 10 10 11 



L B (h;|L ) L B (h- 7 *L Q ) 

Figure 1. Correlations of Sersic parameters of elliptical and SO galaxies. Ellipticals in the Coma and Abell 496 clusters (filled triangles) and Abell 85 (empty 
triangles) measured by Marquez et al. (2000); ellipticals (filled circles) and SOs (open circles) in the Virgo and Fornax clusters measured by D'Onofrio (2001, 
restricted to the good-fitting cases, i.e. where Q < 3 and Xr < 1-5) and dwarf ellipticals in the Virgo cluster (crosses) and in the m vs. L plot measured 
by Binggeli & Jerjen (1998). These symbols are plotted in blue, purple, green, red, and orange in the electronic version of the Journal. The solid curves are 
the adopted fits, performed on log R e vs. log L (eq. 1161 ) and log m vs. log L (eq. 1171 ). and used to estimate ag (eq. 1181 ) and Iq a? s (eq. 1121 ). while the 
dashed curves represe nt the direct fits of log as vs. log L and log Iq a? s vs. log L and the corresponding relation for R e vs. L (from eqs. tm and tEl ). 
The fits use the data of Marqu ez et alj omitting their Abell 85 values, combined with the data of D'Onofrio. The dotted broken line in the m vs. L plot is the 
mean relation from Prugniel & Simien (1997), while the dashed-dotted line is from Graham & Guzman (2003). The luminosities of the D'Onofrio galaxies 
are Ho -independent (see text). 



luminosity relation that disagrees with those of iMarquez et al] and 
iD'Onofriol 

Dis carding the Abell 85 data as well as the dwarf ellipticals 
of lBinggeli & Jerie one is left (circles and filled triangles) 

with fairly strong correlations in Fig. Q whe re, in particular, t he 
shape parameter m correlates with luminosity JCaonet all 1993). 

For simplicity, we assume that elliptical galaxies constitute 
a one-parameter family, based on luminosity (or equivalently, on 
Sersic shape m). Whereas such a one-parameter model of el- 
lipticals is consistent with the Faber- Jackson 119761) relation be- 



tween central velocity dispersion and luminosity, it is obviously 
an oversimplification, in view of the fundamental plane of ellip- 
tical galaxies, where the central velocity dis persion is a function 
of both luminosity and s urface brightness iDressler et aljfl987l 
iDiorgovski & Davisl 19871) . 

We fit the parameter correlations of Fig. Q with 2nd-order 
polynomials in log space (with iterative 3-cr rejection of out liers), 
using the galaxies in Coma and Abell 496 of Marquez et al. (filled 
triangles in Fig.0 and the ellipticals and SOs in Virgo of D'Onofrio 
(circles in Fig.0. Now as and R e are directly related through m 
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R e = b m a s (14) 
b ~ 2m- 1/3 + 0.009876/m , (15) 

where the latter relation is from lPrugniel & Simienl Jl997t) . Simi- 
larly, L and Io a% are directly related through m (eq. 1121 ). There- 
fore, we choose to obtain as and Jo as through the fits of R e vs. 
L, and m vs. L: 



log hro iii fit) 



log m 



(fit) 



0.34 + 0.54 log Lio + 0.25 (log L w ) , (16) 
0.43 + 0.26 logLio -0.044 (logLio) 2 , (17) 
-.(fit) 



where Lui = h'm Lb/ (lO lo L0^, Re"' is measured in kpc, and 
with H = 100 h = 70 fr 7 o kms -1 Mpc" 1 . Then equations ITU 
and I15i lead to 



as 



[6, 



.(fit)^ 



(fit) 



(18) 



The fits of equations 1 16t . 1171 and the corresponding rela- 
tion for as from equation i 1 Si are shown as solid curves in Fig.Q 
whereas the dashed curves show the direct fits to L for as and 
7o 05. For example, for Lb = 10 10 h^XjQ, one has m/^ — 2.7, 
slightl y lower than the value of m = 3.29 inferred from the rela- 
tion of IPrugniel & Simienl < 19971) . The fits can be considered to be 
uncertain to a factor 2 in R e and 1.5 in m and we will propagate 
these uncertainties later in our analysis. Note that equation J 1 61 
produces effective radii that increase at decreasing luminosity for 
low L, which may not be very realistic. Similarly, equation J 171 
produces abnormally low values of m at low luminosity, as the m 
vs L relation shou ld flatten out at low luminosities to accommodate 
the dwarfs of Binggeli & lerien and because values of m < 1 are 
rarely mentioned in the literature if at all. However, inspection of 
Fig. Q indicates that equations <16> and <17t are both reasonably 
accurate for L B > 2 x 10 9 /tf 2 L Q . 



2.2 Scalings of global properties 

We adopt a fiducial luminosity of Lb = L* B (the luminosit 



at the break of the field galaxy luminosity function). Liske et al. 
(2003) have compiled various measurements of the correspond- 
ing absolute magnitude M* and converted to their Bmgc band. 
Their median value is Af MGC = — 19.6 + 5 log h = —20.37 
(for h = 0.7), which is within 0.01 magnitude of t heir co nver- 
sion of the 2dFGRS value of iNorberg et all <2002l) . Liske et alJ 
give the transformation from the Bmgc system to the Landolt 
B system: Bmgc = B - 0.14 5 (B - V). 2 Given th at elliptical 
galaxies have B — V ~ 0.96 iFukugita et alj ll995). we derive 
M B = -20.24, which corresponds to L* B = 1.88 x 10 10 L e , 
again using M® = 5.45 (eq. 1131 ). 

Given the mean luminosity density of the Universe j, and the 
mean mass density of the Universe f2 m p C rit, the mean mass-to- 
light ratio of the Universe is 



3 QtuHq 



(19) 



8nGj 

Converting Liske et al} s j bj = 1.99 x 10 s h L Q Mpc" 3 to the 

B band with B - bj — 0.28 (B - V), initiall y proposed 

by iBlair & Gilmord 1 19821) and inferred from both iLiske et all 



2 The issue of what is the standard B band is confusing, as Lis ke et all 's 
conversion from -Bmgc to ^R.C3 i s quite different from their conversion 
to the Landolt B system. 



and iBlanton et alJ |2003l) . now assum ing (B - V) = 0.94 
jNorberg et all2002l) . with M® = 5.3 jLiske et all) , and convert- 
ing the Hubble const ant, we derive j s / (10 s /170 Lq Mpc -3 ) = 
1.25. Similarly, Blant on et alJ y003) give blue band luminosity 
densities of — 14.98 — 2.5 log h and —15.17 — 2.5 log h magnitude 
per Mpc 3 , when k-corrected to 2 = and z = 0.1, respectively. 
For the local Universe (2 = 0), this translates to jb = 1-04 x 
10 8 h 70 L© Mpc" 3 , again using equation 1 1 31 . Since ILTske et alJ 
do not k-correct to z — and since the median redshift of the 2dF- 
GRS is cl ose to 0. 1, we a dopt a mean k-correction similar to the one 
found bv lBlantonetall -14.98 + 15.17 = 0.19 mag. This then 
leads to a 2dFGRS luminosity density, transformed to the B band 
and k-corrected to z = of jb = 1-05 x 10 8 hro Lw Mpc~ 3 , in 
remarkable agreement with the SDSS value of Blanto net all Tak- 
ing the mean of these 2dFGRS and SDSS luminosity densities, we 
derive a local closure mass-to-light ratio of Yf losnle /hro = 1300 
and, with Q, m — 0.3, we obtain 



T B = 390/iT 



(20) 



So for Lb = L* B — 1.88 x 10 10 h 70 2 Lo, an unbiased universe 
would yield a total mass within the virial radius log hro Mt = 
12.87. 

Two recent statistical analyses of galaxy properties suggest 
that M/Lb has a non-monotonous variation with mass (or lumi- 
nosity), with a minimum value ar ound 100 for luminosities around 
10 10 L jMarinoni & HudsorfeOoilYang et all2003l) . Theoretical 
predictions by Kauffm ann et alj «1999"). based upon semi-analytical 
modelling of galaxy formation on top of cosmological iV-body 
simulations also find a minimum of Tb/Tb = 0.19 at a similar lu- 
minosity, which translates to Tb = 74 hro with equation 1201 . On 
the other hand, the internal kinematics of galaxy clusters are con- 
sistent with th e mass-to-light ra t io of t he Universe given by equa- 
tion (20): e.g. lLokas &Mamonl l2003h derive T s = 351 h 70 for 
the Coma cluster. 

In the following section, we will compare our mass-to-light 
ratios to those expected for old stellar populations. For example, 
iTruiillo. Burkert. & Belli j2004) report a stellar mass-to-blue-light 
ratio T«, b = 7.1 from an analysis of SD SS galaxies using the 
PEGASE jFioc & R gcca-Volmerang el 1997ft stellar population syn- 
thesis code and a Salpeter 1 1955) initial mass function (IMF). How- 
ever, using different stellar population synthesis codes and several 
choices for the IMF, Gerhard et alJl200ll) find values differing by a 
factor of 2 for given galaxies (e.g., T»,s = 4.2 to 8.8 for the nearby 
elliptical, NGC 3379). We performed our own tests of the expected 
values of T„,s given the extinction-corrected colours of nearby 
ellipt ical galaxies, with PEGASE, GALEXV <Bruzual & Charlotl 
2003) and G. Worthey's on-line model 3 . We found a very large 
variety of stellar mass-to-light ratios for given colours, spanning 
from T„,b = 2 to 12, depending on the IMF, its low- and high- 
mass c utoffs, the metallicity and the stellar evolution code. In par- 
ticular. [SalpeteJ 1 s IMF produces high values (as also found by 
iGerhard etaUK and GALEXV produces lower values for a given 
IMF in comparison with PEGASE. In the end, we adopt best values 
of T^b = 8 and 5, and a possible range of a factor 1.5 around the 
geometric mean (T^b — 6.3). 

We also make use of the F aber- Jackson 1 19761) relation. W e 
have compared the calibrations o fide Vaucouleurs & OlsonH 19821) . 
iBender et alii 199(1) . iKormendv et aiT^^TT who make use of data 



3 http : / /astro . wsu ■ edu/wor they/ dial /dial_a_model .html| 
see lWorthe^l 19941 
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from iMcElrovll 19951) andlF orbes & Ponmanl 1 19991 who make use 
of data from Iprugniel & Simienlll996l) . Converting to the same 
value of Ho, we find that the correlation between the veloc- 
ity dispersions avera ged in circula r apertures, a ap and the blue- 
band luminosities of Bender et al] are consistent with the an alo- 
gous c orrelations of de Vaucouleurs & Olson, Kormen dv et alj and 
Forbes & Ponman, but with less scatter, and extends to a large 
enough range of luminosities, with the following mean trend 



171 km s" 



10 10 L 



1/3.25 



(21) 



where Mia = hM„/lO 12 M . 4 In Appendix |X| we derive 
(eq. IA6I ) the concentration parameter for the Nav04 model: 



c = 8.1M, 



-0.11-0.015 logAfi2 



which is similar to the conce ntration paramete rs in equation s Q6J. 

X-ray data analyses JSatoetalJ l2000t IWu & Xud l200Ct 
iLlovd-Davies & Ponma nll2002h based upon equation J2j conclude 
that the total mass of ellipticals, groups and clusters are indeed well 
fit by an NFW model, but their concentration parameters increase 
much faster for decreasing masses. Isato et alj give 



73 Mr 



(27) 



2.3 Distribution of total mass 

We consider here 3 models for the total mass distribution: 1) the 
NFW model with inn er slope — 1; 2) a generalization of the NFW 
model introduced by Jine & Sutd feOOd) with inner slope —3/2 
(JS-1.5, which is slightly different from the formula proposed by 
Moore et al]ll999h : and 3) the convergent model of Nava rro et alj 
( 2004, hereafter Nav04, see eq. Q). The total matter density profile 
can generally be written: 



p(r) 

p{y) 
g{y) 



c*_ 

9(c) 

„— ( 1+1/ )«-3 



(NFW, JS-1.5) , 
exp (-2/iy 1/M ) (Nav04), 



(\n{y + l)-y/{y + l) 



sinb.- 1 ^- y/y/(y + l) 



(NFW) , 



(JS-1.5) , 



(22) 
(23) 

(24) 



k | (2m) 1 - 3 " 7 [3/i,2M2/ 1/M ] (Nav04), 

where —a is the inner slope, c is the concentration parameter 
(eq. (4)), a-h is the radius where the logarithmic slope is equal to 
-2 (NFW, Nav04) or -9/4 (JS-1.5, for which a h /2 is the radius 
where the slope is —2), /i ~ 6 (appendix lAl. and sinh -1 x = 
ln(a; + \/x 2 + 1) for x > 0. In equation 1221 . M v is the mass 
within the virial radius, defined such that the mean total density 
within it is A = 200 times the critical density of the Universe, 
/9 cr it = 3 Hg/ (8 7rG), yielding a virial radius (see iNavarro et alj 
1997) 



2 G M v 
AH* 



1/3 



163 h 1 kpc 



= 206 h 70 1 kpc 




(25) 



for A = 200. Equation 1241 is derived by integrating equa- 
tion 1221 using equation {23j, thus using equation <A2t for the 
Nav04 model. Uing & Suto (2000) measured the concentration pa- 
rameter c (eq. (4)) from their ACDM (with cosmological density 
parameter Q. m — 0.3 and dimensionless cosmological constant 
JIa = 0.7) simulations, which can be fit by the relations 



10.2 Mr, 08 (NFW) 



4.9 M 



12 
■0.13 



(26) 



while for the other two studies, the concentration parameter extrap- 
olated to Mv = 10 12 /i _1 M (Mia = 1) is slightly oyer 100 and 
the slope is —0.51 in both cases. Note that ljing & SutollSato et all 
and Lloyd-Davies & Ponman all define their virial quantities at 
the radius where the mean density is 200 times the critical value, 
while Wu & Xue define their virial radius where the mean den- 
sity is 100 times the critical value, in accordance with the spherical 
top hat infall model for fi m = 0.3, Q a = 0-7. Inte restingly, in 
their cosmological A-body simulations. lBullock et alj 1200 ll) con- 
firm high concentration parameters for the mass distribution within 
subhaloes of haloes (which may correspond to elliptical galaxies 
within clusters), with a steeper dependence of concentration on 
mass: dlnc/dlnM w —0.3. 

The cumulative mass of the dark models used here can all be 
written 

M(r) = M v M(r/a h ) , (28) 

M{y) = (29) 
9(c) 

where g is given in equation J24> . 5 



3 RESULTS 

3.1 Local mass-to-light ratios 

The simplest check that cuspy ACDM models represent the total 
mass distribution of elliptical galaxies, i.e. the gravitational poten- 
tial, is by checking that the mass at all radii is greater than the 
known contribution from stars. 

Given equations J]], 0, J7J, 0, J22}, |23}, and 03, the local 
mass-to-light ratio is 



T(r) 



P(r) _ Pi pjr/a h ) 
Z(r) li £( r /as) 
mT[(3 



p)m] M v 3 _ 

g(c) l v 



xx p - a (l + rix) c 



-3 I II 

exp (x 



(x 1/m ) 



(30) 



(JS-1.5) 



where the second equality is restricted to the NFW and JS— 1 .5 mod- 
els, with x — r/as,Tj = as /ah, and a — 1 (NFW) or 3/2 (JS- 
1.5). 

Fig. |U shows the local M/L profile for the Sersic model with 
luminosity Lb = £* = 1-88 x 10 10 h 70 2 Lq and our three adopted 
mass models with masses within the virial radius of log h^Q 1 M v = 



4 The factor two difference in NFW and JS-1.5 concentration parameters 
arises from the different definition of scale radius is these two models. 

5 Note that the definition of g is the inverse of the definition of g given by 
lLokas & Mamonl 1200 ill for the NFW model. 
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Figure 2. Local mass-to-light ratio profiles for the Sersic distribution of lu- 
minosity L B = L* B = 1.88 X 10 10 h~g L , with m = 3.1 (eq. EH) 
and R e = 3.2h^kpc (eq. EH), and NFW (Ais/ierf cHrvar), JS-1.5 
(dotted cun'es), and Nav04 (ioZtti carves) total mass distributions with 
logh'l M v = 11.87, 12.87, and 13.87 (mass-to-light ratios T B = 
39,390,3900, i.e. 0.1, 1, and 10 times the universal value of eq. l20l ), 
increasing upwards (red, green, and blue in the electronic version of the 
Journal). The open circles represent 0.03 r v , the minimum radius where the 
NFW and JS-1.5 models are accurate in representing the density profiles in 
structures found in cosmological iV-body simulations, while the filled cir- 
cles represent 0.01 r v , which is the analogous radius for the Nav04 model. 
The curves are obtained with equation 1301 using equations EJ. HI. 1141 . 
EH, EH, EU. E3, and EU- The shaded region (orange in the electronic 
version of the journal) shows the mass-to-light ratio of the stellar popula- 
tions of elliptical galaxies (in the range of radii where the Sersic fits to the 
surface brightness profile are believed to be good, see Sec. 12. 21 . 

11.87, 12.87, and 13.87, i.e. for mass-to-light ratios 0.1, 1, and 10 
times the universal value (eq. 1201 ). The shaded region in Fig. 2 
indicates the mass-to-light ratios inferred from stellar population 
synthesis (see Sec. l2.2> . in the range where the surface brightness 
profile is well known (and where the Sersic deprojection is con- 
sidered accurate, see Sec. 12. It . Fig. |2| indicates that, at the lower 
limit of its spatial resolution (0.01 r v ~ 0.6 R e ), the low mass 
(Y_b = 39) Nav04 model produces a lower mass-to-light ratio than 
stellar, reaching values as low as unity, inconsistent with all esti- 
mates of the stellar mass-to-light ratios of old (and red, as observed) 
stellar populations. 

For the models with higher virial cumulative mass-to-light ra- 
tio, the local mass-to-light ratios are consistent with the expected 
stellar mass-to-light ratios in the range where these dark matter 
models are measured. However, simple inwards extrapolations of 
these mass models rapidly lead to sub-stellar mass-to-light ratios. 
For example, for the universal virial mass-to-light ratio (Ts = 
390), the Nav04 model is consistent with the stellar mass-to-light 
ratio at the limit of its resolution, r = 0.01 r v ~ 1.4 R e . But if 
one wishes to ensure that the local mass-to-light ratio at lower radii 
remains higher than the stellar value, one then requires the mass 
model to have a density slope equal to that of a Sersic model. From 
equations ^6j, and |<J. the slope of the deprojected Sersic profile 



Figure 3. Same as Fig.EJfor high concentration parameters obtained from 
X-ray analyses by Sato et al. (2000) for the NFW model and for Lg = 
4.1 X 1O 1O L0 (the median luminosity of the 3 elliptical galaxies in their 
sample), yielding m = 3.7 from equation 1171 and R e = 5.8 h^ kpc 
from equation 1161 . The log masses are now 11.20, 12.20, 13.20 and 14.20 
(respectively corresponding to T = 3.9, 39, 390, and 3900, i.e. 0.01, 0.1, 1, 
and 10 times the universal mass-to-light ratio of equation 1201 , and respec- 
tively black, red, green, and blue in the electronic version of the Journal) 
going downwards at small radii. 



is — p — b/m(r / Re) 1 ^™ . For m = 3.1, one gets a slope of —2.9 at 
1.4 i? e , and of — 1.4 at 0.03 i? e . In comparison, for elliptical galax- 
ies with hM = 1O 12 M0, for which the Nav04 resolution limit is 
r = 0.01 r v = 1.4 R e = 0.08 r_ 2 (eq. (A6|), then according to 
equation the Nav04 mod el has a slope of — 1.3, as also inferred 
from fig. 3 of Na varro et all While the density profile of the Nav04 
model becomes shallower with decreasing radius (as do the density 
profiles of the NFW and JS-1.5 models), one must, on the contrary, 
assume that the profile suddenly steepens at r < R e . The same ar- 
gument applies for the NFW model. Hence, the NFW and Nav04 
models cannot represent the total mass. 

The cuspier nature of the JS-1.5 model makes it less sensitive 
to this criterion. Nevertheless, Fig. [^indicates that a slope steeper 
than —1.5 is required to keep the mass-to-light ratio above the stel- 
lar value. Therefore, unless there is a sharp break in the density 
profiles inside the resolution radius to a slope that is considerably 
steeper than —3/2, any reasonable extrapolation of the mass-to- 
light ratio curves will lead to values lower than stellar. 

Do the high-concentration NFW models found bv lSato et alJ 

(2000) for the total matter distribution also produce abnormally 
lo w mass-to -light ratios? For the three elliptical galaxies studied 
by ISato et all NGC 1399, NGC 3923, and NGC 4636, we find 
from the LEDA database extinction-corrected total blue magni- 
tudes of 10.25, 10.38 and 10.22, respectively, and distance moduli 
from the surface brightness fluctuation study o f|Tonr^UriJ l200ll) . 
corrected by —0.16 magnitude (see Jensen et al. 2003) for th e 
newer Cepheid distance normalization of Freedman e t alJ 1200 lh . 
of 31.34, 31.64 and 30.67, respectively, hence L B = 4.1, 4.8 and 
2.3 x 1O 1O L , respectively. 
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We plot in Fig. [3] the local mass-to-luminosity ratios for 
NFW models with the high concentration parameters of Sato et al. 
(eq. 1271 ) for the median blue luminosity Lb =4.1 x 10 10 Lq. 

The local mass-to-light ratio profiles are increased in the in- 
ner regions, relative to the analogous profiles for the NFW model 
with the low concentration parameters as found in the cosmologi- 
cal JV-body simulations. This is especially the case for low masses, 
which have the most centrally concentrated mass-density profiles. 
Nevertheless, the local mass-to-lig ht ratio found with the high con- 
centration parameter of Sat o et alj is below the stellar value at all 
radii below 0.08 R e , for the very large range of masses considered. 
But here, one could imagine dark matter profiles with somewhat 
steeper inner slopes well within 3% of r v that would produce local 
mass-to-light ratios in excess of the stellar value at all radii where 
the Sersic model is a good fit to the surface brightness profile of 
ellipticals. The same conclusions are reached for the even higher 
concentration param eters ex trapolated to the masses of ellipticals 
from the relations of IWu & Xud ( 20001 and Lloyd-Davies & Pon- 
man. 

To summarise, the low-concentration matter distributions 
found in cosmological N-body simulations produce local mass-to- 
light ratios lower than expected for stellar populations, unless there 
is a sharp steepening in the slopes of the mass density profiles just 
within an effective radius, while high-concentration NFW profiles 
would produce large enough local mass-to-light ratios if their inner 
slopes were slightly steeper. 

These conclusions are unchanged if we vary the effective ra- 
dius or the Sersic shape (rn) by a factor of 2. 

3.2 Velocity dispersions 

Another way to check the compatibility of NFW, JS-1 .5 and Nav04 
potentials with the Sersic luminosity profiles of ellipticals is to 
compute the central stellar line-of-sight velocity dispersion, aver- 
aged within an aperture or a slit. In appendix [B] we re-derive the 
stellar radial and line-of-sight velocity dispersion profiles, assum- 
ing isotropy, and derive for the first time the velocity dispersion 
profiles averaged in circular apertures and thin slits in terms of sin- 
gle quadratures of the tracer density and total mass profile. 6 

Fig. shows the resulting aperture and slit-averaged velocity 
dispersion profiles for Sersic tracers and NFW, JS-1. 5 and Nav04 
mass profiles. The figure clearly shows that, if the total matter is 
represented by an NFW, JS-1. 5, or Nav04 density models, with low 
concentration parameters as found in cosmological simulations, 
then the central aperture and slit velocity dispersions would be 
much smaller than observed. Indeed, the aperture and slit-averaged 
velocity dispersions both obey a v < 113 kms -1 for R < 0.1 R e 
(where the equality is reached for the more favourable JS-1. 5 
model), whereas one expects a v = 208 km s -1 , from the Faber- 
Jackson scaling relation, recalibrated by Bender et al. (eq. 1211 ). 
For the most realistic Nav04 model, the velocity dispersions are 
3.5 times too low for R < 0.1 Re- 
One might worry that this conclusion is reached through heavy 
extrapolation of the dark matter models at small radii. However, it 
is again difficult to imagine a dark matter model that would pro- 
duce aperture or slit velocity dispersions that rise sharply (factor 
of 3 for the Nav04 model from 2 to 0.05 R e ) at increasingly small 

6 One can find in the literatur e the use of tri ple integrals to compute aper- 
ture velocity dispersions (e.g., iBorriello et all2003l) whereas the single in- 
tegrals given in appendix |S]are much simpler to evaluate. 



radii. Given the difference between the mass models of inner slope 
-1 (NFW) and -3/2 (JS-1. 5), Fig.|4]suggests that an inner slope 
< — 2 is required for the mass distribution to recover the high ob- 
served central velocity dispersions. 

Note that the slit-averaged velocity dispersions are slightly 
smaller than the aperture-averaged ones, as is indeed expected 
given that, for the luminosity and mass models considered here, the 
line-of-sight velocity dispersions increase with radius for R < R e , 
and the slit velocity dispersions are less weighted to outer radii than 
are the aperture velocity dispersions. Hence, the slit velocity disper- 
sions are in principle more constraining against a global NFW, JS- 
1.5 or Nav04 potential. However, our quasi-analytical expressions 
for the slit velocity dispersion neglect the effects of seeing and are 
thus of little use at small radii. Moreover, at large radii observers 
subdivide their slit into rectangular bins, whose modelling is be- 
yond the scope of this paper. We therefore, will focus on aperture 
velocity dispersions. 

The high concentration parameters derived from X-ray obser- 
vations produce higher central aperture and slit velocity disper- 
sions, but these velocity dispersions are still well below the pre- 
diction of the Faber- Jackson relation: the NFW model produces slit 
velocity dispersions that are 1.5 times too low for R < 0.1 R e for 
all mass-to-light ratios below 3900 (10 times the universal value of 
eq. EDI)- 

As mentioned in Sec. 12.11 the Sersic shape-luminosity relation 
is probably uncertain by a factor of perhaps 2 (Fig.0. Fig.|5|shows 
the dispersion profiles when the Sersic shape m is taken to be twice 
(top plots) or half (bottom plots) the val ue given in eq uation <17t . 
The higher Sersic parameters f avoured by Bingge l7& Jerienl 1 1 998 ) 
and Prugniel & Simier] ll997l) produce even lower velocity disper- 
sions, and are thus even more constraining against the idea that 
the global potential of elliptical galaxies can be that of the ACDM 
simulations. And if one adopts a Sersic shape parameter as low as 
m = 1.35, as derived from the Marqu ez et alj data for the galaxies 
in Abell 85, one derives central aperture velocity dispersions that 
are almost large enough, but this requires the enormous mass-to- 
light ratio T fl = 3900, i.e. a total mass of over 3 x 10 13 h^M Q , 
which is larger than that of a rich group. For the case of the high 
concentration parameter found in the X-rays, the aperture velocity 
dispersions are still too small, for all reasonable Sersic parameters. 
Moreover, given that the de Vaucouleurs surface brightness profile 
(m = 4) was derived for ellipticals of typically 0.5 to 1 L«, we 
are very doubtful that a Sersic parameter as low as m — 1.27 is 
representative of fairly luminous elliptical galaxies. 

Therefore, the internal kinematics of elliptical galaxies ap- 
pear to be inconsistent with the total mass distribution as found in 
cosmological N-body simulations, even with the steep inner slope 
of —3/2 or the high concentration parameters found by X-ray ob- 
servers. 

3.3 Large concentration mass models as seen in X-rays 

If the NFW model cannot represent the total density profile of el- 
liptical galaxies, how can one explain that several X-ray studies 
converged on a high concentration parameter NFW total density 
profile? 

The answer may be seen in Fig. [6] which shows the stel- 
lar (eqs. (6), (7), and 1121 ) and dark matter (eqs. 1221 . 1231 . 

1241 . 1251 . and 1261 ) density profiles for the median luminosity 
(L b = 4.1fe7n 2 Lo ) of the 3 elliptical galaxies in the sample 
of ISato et alj l200fj) . with effective radius and Sersic shape pa- 
rameter taken from equations 1 1 61 and <17> . respectively, and for 



© 2005 RAS, MNRAS 000, 000-000 



Total mass profiles of ellipticals 9 



R/R 



R / R e 



ra 150 
6 100 

J 5 50 



0.1 1 10 

i i 1 1 ■ ■ ■ | 1 i i 1 1 ■ ■ ■ | 1 — i i 1 1 ■ i ■ | 1 — rm 



0.1 1 10 

T — I I I I ■ 1 1 1 1 I I I ■ 1 1 ■ | 1 — I I ■ 1 1 1 1 1 



~m 350 
300 
i 250 




- 50 



-i i linn 


i i nun 


I i nun 










350 
300 
250 

200 q 

a 

150 P 

ioo !*r 



50 



R/R e 



R / R 



Figure 4. Stellar isotropic velocity dispersions, averaged over circular apertures of projected radius R {top, eq. IB 1 1 1 > or thin slits of half-projected size R 
(bottom, eq. IB 151 ) for Sersic luminosity profiles and CDM dark matter profiles: NFW (dashed curves), JS-1.5 (dotted curves) and Nav04 (solid curves). Left: 
low concentration parameters from cosmological simulations with total luminosity Lb = L* B = 1.88 X 10 10 /i^Lq (corresponding to Sersic parameters 
m = 3.1 and R e = 3.2 h^, 1 kpc). Right: high concentration parameters from X-ray observations (eq. 1271 . NFW model only), with total luminosity 
4.1 X 10 10 h'^LQ (m = 3.7 and R e = 5.8 h~g kpc). Equations j9j, El- IHJ. HH> HH and 113 ms used - The cumulative mass-to-light ratios within 
the virial radius rioo are 39, 390 and 3900 increasing upwards in the left plots and downwards in the right plots (red, green, and blue in the electronic version 
of the Journal). Th e shaded regions are the central observed velocity dispersions expected from the Faber- Jackson ll97d) scaling relation, recalibrated by 
Bende r et al]ll996l eq. 1211 ). The open circles represent 0.03 T2qq, the minimum radius where the NFW and JS-1.5 models are accurate in representing the 
density profiles in structures found in cosmological TV-body simulations, while the filled circles represent 0.01 r2oo> which is the analogous radius for the 
Nav04 model. 



the mass that ISato et all derived within the virial radius (M v = 
3.6 ftfo 1 x 10 12 Mq, leading to T s = 75). Here, we fix the virial 
radius (r v = 316 kpc from eq. 1251 ). and we fit the best NFW 
model to the density profile, within the limits shown as vertical 
lines in Fig. |6| (26 arcsec, set by the PSF of the ASCA telescope, 
and 25 arcmin, which is half the field of view of the G1S instru- 
ment on ASCA). With the median distance modulus of 31.34 (see 
Sect l3.lt . these spatial limits respectively correspon d to 2.3 a nd 
129 /ifo 1 kpc, i.e., 0.0074 and 0.41 virial radii. Since lSato et all fit 
their mass profiles with an NFW mass profile, we minimize the 
square differences in log mass over equally-spaced log radii. 

We find a high concentration parameter: c = 35, although 
not as high as c = 49, inferred from the general concentration vs. 
mass relation of|Sato^t 13 (eq. 1271 ). and only half the value actu- 
ally found bv lSato^talT for these 3 galaxies (which were positive 
outliers on their global distribution of concentration as a function 
of mass). Moreover, using T*,s = 5 instead of 8, the best-fitting 
concentration parameter (with the same quality fit as in Fig.|6j is 



reduced to c = 26, still much higher than c = 9.5 inferred from 
the cosmological simulations of |jin^^iuto] (eq. 1261 ) for NFW 
models with the mass inferred by Sato et all Note that the best fit, 
although adequate, is not superb, because the sum of the stellar 
and dark matter components presents an inflexion point at the ra- 
dius where the two components contribute equally (R ~ 1.1 and 
1.4 R e for T*,s = 5 and 8, respectively), in contrast with the NFW 
model, whose slope steepens continuously with increasing radius. 

Hence, Fig. [5] tells us that fitting an NFW model to the total 
density represented by the sum of the NFW dark matter component 
and the less extended Sersic stellar component can produce very 
high concentration parameters, but the fits are not excellent. 



4 SUMMARY AND DISCUSSION 

Given what we know about the surface brightness profiles of el- 
liptical galaxies, one can simply conclude that the observations of 
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Figure 5. Same as the top plots of Fig.|4] for Sersic shape parameter equal to twice (m = 6.2 and 7.4, upper left and right plots, respectively) or half 
(m = 1.55 and 1.85, lower left and right plots, respectively) of the values derived from eq. 1171 . 



elliptical galaxies cannot allow NFW total density profiles with 
the fairly low concentration parameters as found in cosmological 
TV-body simulations: otherwise one would end up with mass-to- 
light ratios far below stellar within an effective radius, and aper- 
ture and slit velocity dispersions would be far below the predictions 
from the isotropic Jeans equations. Trying a density p rofile with a 
stee per inner slope as a dvocated bv lFukushige & Makinol Il997t) 
and Moore et alTll999l) . using the model of Ijing & SutoThoOO). 
does not help, nor does the recent model by Navarro et al. 1 2004), 
nor probabl y the very similar mass mo del proposed by Stoehr and 
colleagues Istoehr et all2002tlStoehrl2005t) . 

This conclusion of non-conformity of ACDM mass models 
with the Sersic luminosity profile of elliptical galaxies is based 
upon the extrapolation of these mass models to radii smaller than 
the spatial resolution of the simulations. Nevertheless, the only way 
to bring the simulated mass models in conformity with the observa- 
tions is for the slopes of the mass density profiles to sharply break 
around the effective radius from values ~ —1.3 just outside R e to 
steeper values < — 2 inside R e , and it is hard to imagine that this 
sharp break in slope happens at the precise spatial resolution limit 
of present-day cosmological TV-body simulations. This conclusion 
is independent of any assumed value for the stellar mass-to-light 
ratio (as are the curves in Figs.|2|and|3| as well as all of Figs.|4]and 

m 

However, allowing for ten times larger concentration parame- 



ters, as motivated by several X-ray studies, makes the mass-to-light 
ratios and aperture/slit velocity dispersions almost in line with the 
observations, but not quite. 

Assuming no sharp break in the slope of the total mass density 
profile to steeper than —2 inside R e , we are left with an inconsis- 
tency of too low mass-to-light ratios and velocity dispersions of 
the inner regions. This inconsistency can disappear if we add more 
mass in these inner regions. 

The simplest way to add more mass inside an NFW-like mass 
distribution is to postulate that this distribution only applies to the 
dark matter component and that the baryons, essentially stars, dom- 
inate the mass profile in the inner regions, as their Sersic distribu- 
tion has a steep er inner slope than the NFW-like models. In the 
following paper I Mamon & Lokas 2005), we study in more detail 
multi-component models of elliptical galaxies. 

Moreover, during the formation process of ellipticals, the 
NFW-like dark matter should dynamically respond to the dominant 
inner baryons. This process is often studied through the approxi- 
mation of adiabatic contraction iBlumenthal et al. 1986). However, 
the dark matter response to the dominant presence of the stellar 
component will diminish as its local density approaches that of the 
stars. Therefore, the dark matter cannot dominate the stars inside, 
although it could conceivably reach non-negligible densities rela- 
tive to the stellar component. In any event, as long as the inner 
slope of the dark matter density profile in elliptical galaxies is shal- 
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Figure 6. Density profiles: stars (dotted curve, Lb = 4.1 x 10 10 H^Lq, 
m = 3.7, R e = 5.8 h~g kpc, T* t g = 8), NFW-dark matter (dotted 
curve, M v = 3.8 X 10 12 h~^M Q , c = 9.5) and total (solid curve) and 
the NFW fit to the total density (dashed curve), with c = 35. The vertical 
lines delimit the region of the fit (see text). The top abscissa are relative to 
f"200- 

lower than —3/2, it should be difficult to constrain its precise value 
in elliptical galaxies. 

The very different mass distribution for the baryonic (mainly 
stellar) component has important ramifications. It suggests that 
present-day elliptical galaxies have not recently formed through 
rapid collapse, during which violent relaxation should lead to an 
efficient mixing of the baryon and dark matter components, espe- 
cially if their initial distributions were similar (in the linear regime 
of small density fluctuations). 

One alternative is that ellipticals form early by collapse, and 
the dissipative nature of the baryons slowly leads to a segrega- 
tion with the dark matter, which in tur n slowly responds to the 
contracting baryon s by contracting itself. Gned in et alj 120041) and 
|PekeletalJ i2005) find dark matter density profiles scaling as r 
over a very large range of radii in their haloes and elliptical-galaxy- 
like merger remnants of their cosmological and two-galaxy N- 
body simulations, respectively. Although these findings may be in- 
consistent with several observational constraints of elliptical and 
spiral galaxies (perhaps because of insufficient feedback from star 
formation and AGN that pumps energy into the inner baryons, 
which should in turn puff up the inner dark matter distribution), 
they confirm our finding that the total mass density profile cannot 
be as shallow as predicted by the dissipationless ACDM mass mod- 
els. 

The other important alternative is the formation of elliptical 
galaxies (i.e. their morphologies, not their stars, which are thought 
to have formed earlier) through major mergers of spiral galaxies. 
Since spiral galaxies are initially gas-rich and form by dissipational 
collapse, a baryon / dark matter segregation sets in as the baryons 
settle in a disk, while the dark matter extends further out. When spi- 
ral galaxies of comparable mass merge into ellipticals, the baryons, 
which are more tightly bound than the dark matter particles, will 



end up in the inner regions of t he elliptical merger remnant (e.g., 
iDekel. Devor. & Hetzronl2003tlDekel et all2005l) . 

The second conclusion of this work is that considering a two 
component model for elliptical galaxies, summing up the mass 
density profiles of the Sersic stellar component (with a reason- 
able mass-to-light ratio) and a dark matter component as seen in 
the ACDM dissipationless cosmological A-body simulations, one 
finds a total mass density profile that resembles an NFW model, 
with a concentration almost as high as deduced from X-ray obser- 
vations, much higher than in ACDM haloes. Admittedly, the fit is 
not excellent. 

This high concentration NFW total mass distribution is caused 
by the more extended nature of ACDM haloes in comparison with 
the stellar distribution of elliptical galaxies. It would be good if the 
X-ray observers attempt to fit NFW or better Nav04 models to the 
dark matter component of their elliptical galaxies or groups and not 
to the total mass density profile. 

Note that recent modelling by P ratt & Arnau cl l2005l) of their 
XMM-Newton X-ray observations of poor and rich clusters of 
galaxies leads to the conclusion that the total mass density pro- 
file is consistent with NFW models of normal concentration, with 
the same shallow dependence on mass as seen in the cosmological 
A-body simulations. This result is in stark contra st with the strong 
dependence of concentration with mass found by Sato et al. 1 2000) 
from ellipticals to poor clusters to rich clusters. 

One may venture that elliptical galaxies have a more promi- 
nent baryonic component in their inner regions than do groups 
and clusters of galaxies, and will thus, by superposition of Sersic 
and ACDM components, appear more concentrated when fitted 
with single NFW models, than expected from the general trend of 
groups or clusters. However, ISato et al. find a continuous power- 
law trend for the mass dependence of their concentration param- 
eters. In other words, Sato et al. find moderately high concentra- 
tion parameters for groups of galaxies, as do Wu & Xue 1 2000) and 
iKhosroshahi et al] 120041) . suggesting an important baryonic contri- 
bution in their inner regions, albeit less dominant than in elliptical 
galaxies. 

Interestingly, our best NFW model for the dark component 
of the 2-c omponent model we try for the galaxies observed by 
ISato et all has a conce ntration, c ~ 26 to 35, that lies in between 
the very high value of ISato et all c ~ 49, and the low value ex- 
pected from cosmological A-body simulations for objects of that 
virial mass, c ~ 9.5. Additional measurements of concentration 
from X-ray observations of ellipticals should greatly clarify this is- 
sue. 

1 1 the following paper l Mamon & Lokas 120051) . we investigate 
in more detail a 4-component model, including stars, dark matter, 
hot gas and a central black hole, and ask whether the observations 
of ellipticals require cuspy cores, and to which accuracy one can 
measure the total mass within the virial radius through kinematic 
observations. 
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APPENDIX A: CONCENTRATION, CENTRAL DENSITY AND OUTSIDE MASS VERSUS VIRIAL MASS FOR THE 
NAVARRO ET AL. (2004) MODEL 

iNavarro et alj 120041) have recently shown that the density profiles of haloes in cosmological iV-body simulations begun with a ACDM 
spectrum can be fit to high precision to 



p(r) = p-2 exp(2/tt) exp 



-2fi 



r~2 



(Al) 



where our p is the inverse of their a and where p2 is the local mass density at the radius r_2 where the logarithmic slope of the density is 
-2. 

The enclosed mass of the density profile of equation <A1> is 



M(r) = 4tt p-2 rt 2 2" 3m p 1 ' 3 ^ exp(2 p) 7 



3p, 2p 



r-2 



(A2) 



where y(a, x) = J* t a_1 exp(— t) dt is the incomplete gamma function. INavarro et alJ i2004h provide a table of values of 1/p, which for 
giant galaxy mass objects yield a geometric mean and median of (p) = 5.85 and 6.16, respectively (and roughly the same for haloes with 
dwarf galaxy or galaxy cluster masses). We therefore adopt p = 6. 

The reader may note a resemblance to the enclosed luminosity of the Sersic m odel given in eq uations J 101 and i 1 1 i . which may be 
related to the resemblance of the projected NFW model to the Sersic profile noted by Lokas & Mamon 1 2001 ), and the value of p is of the 
rough order of the Sersic m's (see Fig.0. 

Expressing the mean density at the virial radius r v = cr-2, one obtains 



Bzl = 1* (2p) 3 ^ exp(-2 M ) , C " r/ - 



= /i(c) 



Now. INavarro et all provide a table with p-2/pcilt an d r_2. We find that their data can be fitted by 



P-2 
Pcrit 



dex 



4.37 - 0.30 log 



hr. 



1 kpc 



0.10 loe 



hr- 



1 kpc 



h 



hr~ 



1 kpc 



(A3) 



(A4) 



We solve this 2nd order polynomial for /ir_2/(lkpc) = f 2 {fi-il 'Pcrit)- Equations <A3t and <A4> can be combined into an implicit 
equation for the concentration parameter c: 



c/ 2 1 [/1(c)] = r v 



(A5) 



Given that the virial radius is a function of the mass M v within the virial radius (eq. 1251 ). then for a given M v , one can solve equa- 
tion <A5> for the concentration parameter. For p = 6, we obtain 



1M 



-0.11-0.015 log A/12 



(A6) 



where M12 = /iM„/10 12 M s . 

Fig. lAll displavs a check of equation <A6t : for every M v , we obtain c from equation <A6t . p-2/pcrit from equation <A3t . but also r v 
from equation I25i and r_2 = r v /c. The agreement is excellent and much better than if we had instead directly obtained f^ 1 from a fit of 
/ir_ 2 /(l kpc) as a function of p-2/ Pcrit- 

The central density of the Nav04 model is 



pa = p(0) = p-2 exp(2p) 

Using equation <A3t . equation <A7t can be expressed in terms of the critical density 

Pcrit 3 K ^' 7(3m,2 M c 1 /m) 
and with equation JA6> . this yields 

Po 



log ~ 9.0 - 0.25 log M12 - 0.030 log^ M 12 . 

\ Pcrit J 

Thus the central density decreases with increasing mass. 

The mass of the Nav04 model converges at infinity, and the mass at the virial radius satisfies 



(A7) 



(A8) 



(A9) 



M v _ y(3p,2p c 1 ^) 



0.44 M, 



-0.068-0.014 log Mi 2 -0.0012 log^ M 12 



(A10) 



Moo T(3p) 

so that for giant galaxies a little over half the mass is beyond the virial radius, while for rich clusters more than 3/4 of the mass is beyond r v . 
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Figure Al. Normalised density at scale radius versus scale radius for the Navarro et al. (2004) model. Circles: measurements by Navarro et al.; curve: 
predictions from our concentration parameter / mass approximation (eq. IA6I ^ also using equations 1251 and IA3I : numbers: log hjQ M v . 



APPENDIX B: LINE-OF-SIGHT, APERTURE AND SLIT VELOCITY DISPERSIONS OF ISOTROPIC SYSTEMS 
Bl Radial velocity dispersion 

The Jeans equation is 

« + 2M^ = -^)^M, (Bl) 
dr r r 2 

where the anisotropy parameter is fi — 1 — a 2 /of, with a t = erg — er^, the ID tangential velocity dispersion, so that [3 = corresponds to 
isotropy, j3 = 1 is fully radial anisotropy, and j3 — > — oo is fully tangential anisotropy. 
For isotropic orbits, the Jeans equation IB II trivially leads to 



G 



£M 



(B2) 



B2 Line-of -sight velocity dispersion 

Projecting along the line-of-sight, one trivially finds that the line-of-sight velocity dispersion is 

lol r dr 



I(R)crt os (R)=2 



Inserting equation IB2> into equation <B3> . and inverting the order of integration, one easily finds (Prugniel & Simien 1997): 



I(R)af OB (R) = 2G 



Vr 2 -R 2 



£(r) M(r) dr 



(B3) 



(B4) 



Writing M(r) = M v M(r/a$), x = r/as, x v = r v /a$, rj = as/an, the deprojected luminosity L3 of the Sersic profile is given in 
equation 11 W . the isotropic velocity dispersion, derived from equation <B4> . using equations Q, |6|, 0, J8j, and I28t becomes 



O-los(-R) 



T(2m) 
r[(3-p)m] ^ 



V 2 exp ( X 



l/m\ 



Vx 2 -X 2 



£(x) M(x)dx , 



(B5) 



where V 2 = GM v /r v is the squared circular velocity at the virial radius and again X = R/as- 



B3 Aperture velocity dispersion 

The aperture velocity dispersion satisfies 
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(B6) 



where L 2 (R) = J Q 2nSI(S) dS is the luminosity within projected radius R. Inserting Iaf OB from equation <B4t into equation <B6> . again 
inverting the order of integration, we obtain the isotropic solution for the aperture velocity dispersion: 



4ttG 
MR) 



L 2 

■ Otot 



rl{r)M(r)dr 



4nG 



(r 2_ jR 2 ) 3/2 



e(r)M(r)dr 



£{r)M(r)dr , 



L 2 {R) tm 3L 2 (R) J R 

which, in the limit R — ► oo, converges to the (one-dimensional) isotropic velocity dispersion, averaged over the entire galaxy: 

AtyG 



ftot = 



3L 



rl{r)M(r)dr . 



(B7) 
(B8) 

(B9) 



Using the formula for the projected luminosity of the Sersic profile iGraham & Collessil997tlBinggeli & Jerierll998tlLima Neto et alJ 



L 2 (R) — 2nm Io a s 7 



R \ 1/71 



2m, (*) 
\a s J 



(BIO) 



the isotropic aperture velocity dispersion at radius R = X a$ (eq. IB8I ), together with equations |6|, Q, j8j, <10t . <28> , and <B 101 yields 



2 (m r(2m)/r[(3-p) m] r v a 



a;^(a:)M(a:)da; - 



■ M(x) da; 



(Bll) 



For the numerical integration of equation IB 1 li . the first integral can be evaluated by integrating along In x in the range [—20, In a; max ] , 
where x max = 50 m , for which the exponential term in t(x) is extremely small. The other integrals can be numerically evaluated by 
integrating \n(x/X) in the range [0, ln(a; max /X)]. In this paper, we evaluate numerically all velocity dispersions using Mathematica. 



B4 Slit-averaged velocity dispersion 

The velocity dispersion averaged over a thin slit of width R is 

2 J R i(s)*L(s)ds 



(B12) 



/o I(S)dS 

With J(R) = J* I(S) dS, equation jB12i becomes 

J W°™W = 2 ["dS r 

G Jo Js r 

pR roc j poo / Z?\ 

= I / i(r)M(r)dr + R / £{r)M(r) ^Jr 2 -R 2 + / £(r) M(r) sin -1 ( - ) dr , (B13) 
2 Jo ./« r Jr ^ r ' 

where the second equality in equation IB 13i was obtained after inverting the order of integration in the first equality. With equation Q, one 
has 



J(R) = mloas'f 



(-) 

\a s / 



(B14) 



The isotropic slit velocity dispersion (eq. IB 131 ) becomes, with equations |6), {7), 1281 . and <B14> . 



fslit(-R) = 



r(2m)/r[(3-p)m] r v T/2 



2m7(m,X 1 /'") as 



v i> 



l(x)M{x) dx + X 



^x 2 -X 2 



£(x)M(x)dx+ I sin -1 ( — ) ?(a;) M(z) da; 



'0 

where X = R/as- The last two integrals are easily evaluated numerically with the substitution x = Xj sin it. 



(B15) 
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